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Abstract. We find a gravitational wave solution to the linearized version of 
quadratic gravity by adding successive perturbations to the Einstein's linearized 
field equations. We show that only the Ricci squared quadratic invariant 
contributes to give a difli'erent solution of those found in Einstein's general 
relativity. The perturbative solution is written as a power series in the fS 
parameter, the coefficient of the Ricci squared term in the quadratic gravitational 
action. We also show that, for monochromatic waves of a given angular frequency 
u!, the perturbative solution can be summed out to give an exact solution to 
linearized version of quadratic gravity, for < ui < c/ This result 

may lead to implications to the predictions for gravitational wave backgrounds of 
cosmological origin. 
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1. Introduction 



Gravitational Lagrangians which include quadratic curvature terms, have been 
proposed at classical level as an extension of Einstein's theory More recently, 
quadratic Lagrangians have been used to yield renormalizable theories of gravity 
coupled to matter [2]. Higher derivative theories also arise as the low energy limit 
of string theories |31 ^ . 

For the sake of definiteness we will deal with the quadratic (fourth-order) gravity 
generated by the action 



^= TTTT^ / d^xV^ {R + aR^ + m,.uR'"' + IQttGC^}. (1) 
iDTrG J 

Here, R is the Ricci scalar, i?^^ is the Ricci tensor, G is the Newton's gravitational 
constant, and a, (3 the coupling parameters of the quadratic interactions. We use 

throughout this work units in which h = c = 1, metric signature (H ), Riemann 

tensor i?"/3^,. = -r^^.^ H and Ricci tensor Rfj_„ = g^^Ra^xup- 

Let us consider the action S. For a = /3 = 0, 5 is exactly reduced to the Einstein- 
Hilbert action. Then, we can expect that for a and j3 sufficiently small the quadratic 
gravity can be approximately described by Einstein's General Relativity second-order 
equations with small corrections due to quadratic curvature couplings. In this case 
approximate solutions to quadratic gravity can be constructed by adding successive 
perturbations to a solution of the Einstein's gravity. This perturbative approach 
is develop from the concept of regular reduction of a system of partial differential 
equations and is used by Lousto et al. CI IE] to obtain perturbative solutions 
to quadratic gravity for some non-radiative astrophysical scenarios. In the present 
work instead, we use a perturbative approach to build radiative solutions to quadratic 
gravity. 

In section |21 we give a description of the scheme on how the perturbative field 
equations are obtained. In section |3 we use the perturbative method to build a 
solution to the linearized quadratic gravity. In sectional we discuss the implications 
of the solution obtained to astrophysics and add some considerations concerning the 
detectability of the computed effect. 

2. The general perturbative scheme 

We give in this section, a description of the general scheme on how the perturbative 
equations are obtained. We then show how a solution to quadratic gravity, as a power 
series in a and (3 parameters, can be built. 

Let us now write the general field equations to the fourth-order gravity. Variations 
of the action S, with respect to the metric g^i, leads to the following field equations: 

G^. + aH^, + 131^, = -SttGT^,, (2) 
where Gp,^ is the Einstein's tensor given by 

Gi_ii, ~ Ri_ii, — —gi_ivR- (3) 
The terms in the fourth-order derivatives of the space-time metric are given by 
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Hf,, = 2R^, - 2g^,UR + ^^g^.R' - 2RR^, (4) 

and 

= R,tiiy — ^dt^"^-^ ^ ORf^i, — 2Rfj"Rav + ■^gfivRa/sR"'^ ■ (5) 

Let us now write the perturbative equations for the system that is, the 
equations which give as solutions the contribution to the space-time metric of a given 
order in a and (3 parameters. 

In the zeroth-order (a — (3 = 0), the perturbative equation is by construction the 
ordinary Einstein's gravity equation, which is written as 

GW = -87rGr(«), (6) 

where 

G[.°j = <)-i.g,.i?(°). (7) 

The solution of these equations is denoted by g^jS and is independent of a and (3 
parameters. From the above relations we obtain 

= 87rGr("' (8) 

and 

<^ = - lv,.T("^) . (9) 

At the Ith-order we have the equations 

GW = -87rGr^i) - - /^/^ , (10) 

where H^^J and /^^^ can be written in terms of the energy-momentum tensor from 
and (0). The solution of these equations at first order in the coupling parameters a 
and 13 is denoted by g^^} . 

By repeating this procedure to higher order we obtain at nth- order the following 
of perturbative equations: 

G("j = -87rGr(i;) - aij'r " /?^^r ■ (H) 
These equations show that at each order of perturbation scheme we have the 
Einstein's General Relativity equations with a new effective energy-momentum source. 
Thus, in this perturbative scheme, the quadratic curvature couplings are left out of 
the dynamics, which is given by Einstein's gravity equation. The higher order terms 
enter only as components of the source in the perturbative field equations. 

The approximate solution to the quadratic gravity can be constructed by adding 
the contributions to the space-time metric order by order, that is. 



a - o(0) I I -I- 



(12) 
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3. The perturbative scheme in linear approximation 

The ordinary theory of gravitational emission and propagation is developed in the 
radiation zone, where the variations of the gravitational field can be treated linearly. 
Then, we apply the scheme delineated above to the linearized quadratic gravity. To 
deal with the linearized theory we expand the metric g^i, as follows 

5^1/ = Vtii^ + Khf,^, (13) 

where k <C 1, and 7/^,^ is the Minkowski space-time metric. Within this approximation, 
we can neglect the nonlinear curvature terms present in Q and Thus, the 

linearized version of the eq. Q can be written as 

Glt) = -«.r,.-2^4^)-2^/(^), (14) 
where k = VIOttG, G^il) is the linearized Einstein tensor defined by 

G(,t^ ^ Uh^, - + h^, ~ r^^,{Uh - K^^^p) (15) 

and the terms which involve fourth-order derivatives of the metric reduce after 
linearization to 



(L) 



2Dh^, - - 277^,(n2/i - Dh'^^^^is) (16) 



and 



K ^ z (17) 

Let us adopt the Hilbert gauge, namely 

h^u" = \k. (18) 
and define the trace reverse of h^^ by 

if^^ = /i^i. - ^?7pi//i- (19) 

In our approximation we deal with metric independent sources, therefore the 
energy-momentum tensor is independent of a and (3 parameters. Thus, 

T[a)^T^. and ^ ^ 

7 20 
T^;^' =0 for n > 1. 

Now, we write the perturbcitive equations expanding the field variable ^^u-, 
following the scheme delineated in the previous section. 
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The zeroth- order perturbative equation is the linearized Einstein's gravity 
equation 

□^w = (21) 

whose sohitions in the radiation zone are given by the usual retarded gravitational 
potential 



(0) ^ _ 

47rr 



^y'dVT^,(t-r,x'), (22) 

where r denotes the distance from the source to the observer which is very large as 
compared to the dimensions of the source. 

These are just the ordinary Einstein's General Relativity solutions from which we 
know that they can be made transverse, traceless and have only spatial components. 

By virtue of H2U|I . and the traceless nature of ip^j^j^ the 1-th order perturbative 
equations are independent of a and reduce to 

□^w = -p^nr^,. (23) 

These equations can be solved by applying the usual Green's function method. 
Cancelling the surface integrals we obtain 

which can also be written in the form 

^11^-/3^^2^- (25) 



dt^ 
tion 

The approximate solution up to Ith-order is given by 



Since (/j^]'' is proportional to (/j^^^ the former is also transverse and traceless. 



Since (/^'-^^ = 0, h^^i = </'[p. Thus we can compute the scalar amplitude, 



|_ + /i^ , for the first quadratic curvature order solution via equation l|26(l . 
At Ith-order in /3 we have 

= h(") + /? ""^"y^^ , (27) 

where 



is the scalar amplitude given by Einstein's General Relativity and. 
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For monochromatic signals of angular frequency w we get the simple result 

h.^'^=(i+^y.^'\ (28) 

whenever |(/3w^)/c^[< 1, with c being the velocity of Hght. The above result shows 
the amplitude of a single monochromatic gravitational wave in quadratic gravity to 
the first approximation in the higher curvature terms. 

3.1. An exact solution to linearized quadratic gravity 

By adding the contributions of higher-order perturbative terms, the solution at the 

n-th order can be written as 

^ki =[^ + ^ + ^ + ^ + --- + ^;r)^ki- (29) 

The series in (/3a;^)/c^ is a geometric one. This series can be summed to give the 
exact result (1 — whenever |(/?w^)/c^|< 1. Therefore, an exact solution for 

a monochromatic wave in quadratic gravity in the linear weak field approximation is 

given by 

where /ij^"^ is the Einstein's General Relativity solution. 

The scalar amplitude hs for this solution is related to Einstein's gravity scalar 
amplitude h^^^ by the same factor, i. e., 

hs=(l-^)\^'^. (31) 

We observe that this result is valid only for |(/3w^)/c^|< 1 or, for waves with frequency 
0<w<c/|/?|V2. 

The result derived above shows a correction to the amplitude of linear oscillations 
of the space-time metric due to the presence of a Ricci-squared term in gravitational 
Lagrangian. This result is remarkable simple and predicts a frequency dependent shift 
in the wave amplitude. We point out that, in the context of higher-order gravity, this is 
an interesting result since it enable us to know the effect of quadratic curvature terms 
in linear oscillations of space-time metric by avoiding the negative energy problems 
that can arise due to the effect of a massive ghost, associated to the Ricci squared 
term, on the linearized fourth-order field equations. Since the original fourth order field 
equations are reduced, by the perturbative scheme, to the linearized Einstein's gravity 
equations in presence of an effective source, the energy density and the energy flux 
of gravitational waves can be computed from the same expressions for this quantities 
derived in Einstein's linearized gravity. 
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4. Possible implications to astrophysics and cosmology and perspectives 
to detectability of the amplitude correction 

Let us now discuss how the result derived in this paper can affect the predictions for 
the gravitational radiation generated by astrophysical and cosmological sources and 
discuss a possibility to detect the predicted quadratic curvature amplitude correction. 

Firstly, we need an order of magnitude estimative for the amplitude shift. To first 
order in /3 the shift is given by the product of Einstein's linearized wave amplitude 
by the factor (3(J^ /c^. Observational constraints on /3 parameter derived from the 
bending of light by a gravitational field in the fourth-order gravity indicates an upper 
bound of < 10~'*cm^ [S]. The same bound is obtained form a translation of the 
results of the sub-millimeter tests of the inverse square law of gravitational attraction 
among massive bodies [TSj- Due to the small value of (3 upper bound, the effect of the 
Ricci squared correction to the linearized wave amplitude will appear only for high 
frequency gravitational waves. Even for inspiraling compact binaries of neutron star or 
stellar black hole pairs, which experiments an intense growth of the frequency when the 
system approach the coalescence, the frequency of the gravitational oscillations cannot 
exceed lO^Hz. In this case, the shift in the wave amplitude will be < lO^^^ft.'^*'-', where 
/i'^"^ is the Einstein's linearized wave amplitude. This shift is clearly non observable. 

For the amplitude shift to be observable, taking the maximal value acceptable for 
/3, the frequency of gravitational waves must lie in the range of 10^°— 10^^ Hz. Compact 
binaries formed by pairs of mini black holes having masses M <Ax lO^^M© are able 
to emit gravitational waves with frequencies above 10^" Hz. As pointed out in a 
significant fraction of the dark halo of our galaxy can consist of a large number (~ 10^^) 
of sub-lunar mass primordial black holes with mass lO^^^A/© M < lO^^M©. 
Therefore, it is natural to expect that most of these objects form binaries which 
can emmit very-high-frequency gravitational radiation before coalescing. Then, we 
can expect a background of gravitational radiation associated with these objects in 
the frequency band where the quadratic curvature effects might be relevant. 

On the other hand, gravitational waves with frequencies in the GHz band 
are predicted by several models for stochastic backgrounds of gravitational waves 
of cosmological origin derived from inflation, string cosmology and cosmic strings 
scenarios |12l 113) . Then, we can expect that the result derived in the present work 
can affect the predictions for cosmological backgrounds of gravitational radiation, in 
the high (GHz) frequency spectral band. An indication of the effect of quadratic 
correction on the spectrum of a stochastic background, which is assumed to be 
isotropic, stationary and unpolarized can be easily computed. The energy density 
of the gravitational waves is given by 

^2 . . 

where hki is the time derivative of the transverse traceless gauge filed and < > means 
spatial average over several wavelengths. If we insert the quadratic correction to 
gravitational waves in the equation H32|l and compute the ensemble average for a 
stochastic background we get the result 

n,uf) ^ = f 1 + 8-^^) ^f2if), (33) 

Pcd{logf) V ) ^ 

where / is the wave frequency (uj = 27r/), flf2if) is the "spectral function" predicted 
by Einstein's theory and pc is the critical energy density required to close the universe. 
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We note that the ilg^ wiU differ from il^J only in the (10^° — 10^^) Hz frequency range. 
Such very-high-frequency gravitational waves reflect the behaviour of the very early 
universe jj^]- We note that the expression H33|) must be taken with care since the 
computation of pgw may be different for a given cosmological models. Thus, a deeper 
inspection of how the cosmological background is affected by the quadratic curvature 
corrections must be carried out. 

Once very-high-frequency (10^° — 10^^ Hz) gravitational radiation could be 
detected, this detection could be used to set up new observational constraints on the 
P parameter. The gravitational wave detectors such as LIGO, VIRGO and the mass 
resonant detectors are not adequate to observe such high frequency waves. The effect of 
quadratic curvature gravity on the linearized fluctuations of space-time metric will only 
be observable when a new generation of detectors with sufficient sensitivity to measure 
amplitudes of gravitational waves with frequencies in the range of (10 — 10^) GHz are 
available. 

A gravitational wave detector which explores the effects on electromagnetic 
fields of the passage of gravitational waves has been proposed by Cruise jTJI. The 
proposed device can be designed to observe gravitational radiation with frequencies 
in the MHz and GHz range. Although the sensitivity of the Cruise detector does 
not reach the required sensitivity to detect gravitational wave signals, improvements 
in the technology might lead to more sensitive detectors in the near future. This 
detectors will be able to observe the gravitational radiation generated by cosmological 
sources at the very early moments in the universe. Therefore, the quadratic curvature 
corrections to gravitational waves might be detectable by a sufficiently sensitive 
generation of electromagnetic detectors. These detectors might also be able to observe 
the background radiation produced by a halo population of binary systems composed 
of sub-lunar mass primordial black holes, whereas these objects may actually exists. 
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